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We compute the total cross section as well as angular and energy distributions for process e''"e“ —> 

777 with both unpolarized and polarized beams in the framework of noncommutative quantum 
electrodynamics (NCQED). The calculation is performed in the center of mass of colliding electron 
and positron and is evaluated for energies and integrated luminosities appropriate to future linear 
colliders. We find that by using unpolarized beams it is possible to probe the Lorentz symmetry 
violating azimuthal dependence of the cross section. Furthermore, with polarized beams the left- 
right asymmetry of the CP violating NCQED amplitudes can be used to obtain bounds on the 
noncommutative scale Anc which exceed 1.0 TeV. 


1. INTRODUCTION 


The idea of formulating field theories on noncommutative spaces goes back some time P|. Interest has been 
revived recently with the realization that noncommutative quantum held theories emerge in the low energy limit 
of string theories iSJl. This has led to numerous investigations of the phenomenological implications of 
noncommutative QED |^|^. 

In noncommutative geometries, the coordinates obey the commutation relations 

( 1 ) 

where = —9'^^. The extension of quantum held theories from ordinary space-time to noncommutative 
space-time is achieved replacing the ordinary products with Moyal * products, dehned by 

{f*9){x) = exp{kiS^''da^^dy^)f{x)g{y)\^^y . ( 2 ) 

Here, in order to ensure the S matrix unitarity, we assume that = 0, i = I, 2, 3, . 

In the following, we study the effect of noncommutative geometry on the process e^e“ ^ 77T using non¬ 
commutative quantum electrodynamics, NCQED. NCQED, dehned and described for example in [ 3 , has as its 
Lagrangian 

C = 'tjj-k{p-m)tp - ^Tr + Cgauge + Cghost , (3) 

where Cgauge and Cghost denote the gauge hxing and ghost terms. The corresponding Feynman rules for phe¬ 
nomenological calculations can be derived from Eq. © m Since scale at which noncommutative effects are likely 
to occur is large, we focus on the energy scales typically associated with future linear colliders. Our calculations 
are performed in the center of mass of the colliding electron and positron. In the next section, we outline the 
calculation of the squared amplitudes. This is followed by a summary of the cross section computations and a 
discussion of the results. Details of the calculation are presented in the Appendices. 
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2. NCQED AMPLITUDES 


Typical Feynman graphs contributing to the e+e“ ^ 777 at leading order are shown in Fig.^ Diagrams in 
which the fermion line is connected to the final photons by a single photon propagator vanish. The complete 
set is obtained by permuting the photons, which gives a total of twelve diagrams. We have calculated specific 
helicity amplitudes and it is therefore unnecessary to include ghost contributions The amplitudes were 
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FIG. 1: Diagram (a) represents the contributions from abelian interaction terms and diagrams (b) and (c) represent 
contributions from the non-abelian interaction terms. 

calculated using the Feynman rules in Ref. 0 ■ We computed the helicity amplitudes with the aid of the symbolic 
manipulation program FORM and simplified the results using Mathematica. 

To simplify the presentation of the results, we introduce kinematical variables 


a = k 2 -k^b 

= k3-kic = ki-k2 

(4) 

II 

= pi-k2r=pi-k3 

(5) 

s = p 2 -kit 

= P2-k2U = P2-k3 

( 6 ) 

{ki-9-k2/2) 

w = sin {ki-9-k2/2) 

(7) 

S 

= {Pl +P2f , 

( 8 ) 


where p-9-q = These variables are related as 

S/2 = s + t + u = p + q + r 
a = q+r—s = t+u—p 

b = pr — t = s + u — q ^ ' 

c = p + q — u = s + t — r. 

In the center of mass system we also have 

ki-9-k2 = k2-9-k^ = k'i-O-ki . ( 10 ) 


While there are only five independent variables, it is convenient for displaying the results and examining sym¬ 
metries to retain six. If the helicities are labelled A, A; Ai, A 2 , A 3 , the square of helicity amplitude -|-, —; -|-, -|-, — 
is 








3S o -I-s^) b[q‘^+t‘^) c(r^-l-M^) 
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( 11 ) 


The remaining squared amplitudes are given in Appendix R71 

With the aid of Eqs. (141511)171 and llOII . it is easy to check that Eq. d satisfies Bose symmetry. The last term 
in Eq. d changes sign under the exchange of pi and p 2 , which is a reflection of the lack of charge conjugation 
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symmetry. The parity transformation p —p, X —> —A gives A1 a,A;Ai,A2 ,A3 = as can be 

seen using ApoendixO As a consequence, CP is violated, but CPT is preserved The CP violating terms 
cancel in the sum over fermion helicities. 

Summing over photon helicities and averaging over electron helicities gives 

gps + s^) + qt {q^ +t‘^) +ru (r^ + 
pqr stu 

'^-4^2 / a(p^ + s^) ^ b{q'^+t^) ^ c {r"^ + u^) pg + [p + s) {q + t) + st 
\ be ac ab c 

pr + {p + s) {r + u) + su qr + {q + t) {r + u) + tu 3 S 
b a “ 

We have checked that this result satisfies the Ward identity = 0 for each photon and that it reduces to 

the standard QED result 0 if = 0. 
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3. CROSS SECTION RESULTS 


3.1. Unpolarized Cross Section 


The details of obtaining the cross section from |Ale+e-^ 777 pj Eq- OUi given in Appendix IXI The result 
consists of the pure QED cross section, ■which contains an infrared divergence that must be regularized, and an 
infrared finite NCQED correction. To check the validity of our helicity amplitudes, we recalculated the QED 
spin averaged total cross section by retaining all terms linear in the electron mass squared, in the numerators 
of Eq. 11211 . Our result then agrees with that of Berends and Kleiss ^3; namely 

'^2-1-^777 = ^ ((logp - 1)^ 0-Oguj - 1) + 3) , (13) 


with 


P = 


S 


U! = 


JT' 2 
'^min. 


(14) 


By assuming the non-commutative factor in Eq. 03 is small, we have = sin^p = where ip is 

p = \ki-9-k2 = \{kiX k-2^-9. (15) 

It is then possible to express the non-commutative effects in terms of A, the angle between pi and 9 and 
2 ; = S/Aj^Q. The NCQED contribution to the spin averaged cross section has the form 


a 


NCQED 

■ e+e~ —>777 


r 


S 


2231 TT^ 


/ TT^ 


~ + sin A -I- 2^(3) ~ 


148957 71ogp 


720 120 2 

If we average over A, the expression for the total cross section becomes 

'-3 r _ / 65219 TT^ 3 


34560 


960 


a 


(Jg+ 


g-rg- —>777 


2(logp-l) (logw - 1)-I-6-b z 


_ro3)-I^^ 

V69120 480 4^^ 1920 J 


(16) 


(17) 


With or without the average over A, the effect of non-commutativity on spin averaged total cross section is 
relatively small since it depends on and has a logp rather than a log^ p dependence. 

To determine if NCQED can be tested using the spin averaged three photon process, we examined the depen¬ 
dence of the cross section on the azimuthal angle (f) of one of the photons. This is a pure NCQED effect since 
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the QED cross section has no such (/(-dependence. Before doing this, it is important to determine the range of 
validity of the approximation sin^ ip k . From Appendix^ ip'^ has the form 

sin^ a sin^ (3 , (18) 

54 

where Vi = Ei/E and a and (3 vary between 0 and tt. Thus, the approximation is always good if z < 8 . Since all 
the other factors vary between 0 and 1 , and we integrate over some or all of them in calculating the distributions, 
a limit on z based on the the average value {<p^) can be useful. This limit is 


z < 




< 35. 


sin^ asin^ (3) 


(19) 


In the following, we choose ^fS = 0.5 TeV, 1.0 TeV, 5.0 TeV, which, for Ajvc = 1.0 TeV, corresponds to z = 
0.25,1.0,25.0. 

As in the e+e" ^ 77 case, the ^-dependence of the spin averaged cross section is proportional to z^. With 
no cut imposed on the polar angle 6, we have 


da 

dfj) 


z^a^ [2231 TT^ 5^,^, . 2 , /86141 

[l440 V 69120 ~ 


1963 TT^ 71ogp 
KM) ^ ~ 1440 




IItt^ 

1440 


71ogjO 

1152 


C(3) 


( 20 ) 


The effect of this characteristic (/(-dependence is illustrated in the left panel of Fida The signature of 37 (j)- 




FIG. 2: (color online) In the left panel, the solid line is number of events as a function of (/> for the case y/S = 5 TeV, 
Ajvc = 1 TeV, A = rr/d, C — 500 fb“^ and no cut on cos9. The dashed line is the uniform background from QED with no 
cos 6 cut. The right panel shows the effect of imposing the additional cut | cos 6 | < 0.9. 


dependence, a unique feature of NCQED, can be further enhanced by imposing a cut on polar angle 9. This has 
a rather large effect on the NCQED signal since the QED contribution to the e+e” ^ 777 cross section is very 
sharply peaked in the forward a backward directions. The result of the cut | cos0| < 0.9 is shown in the right 
panel of Fig|21 

We also checked the NCQED corrections to the QED energy and polar angle distributions of one of the photons. 
While there are some differences in the shapes of the NCQED distributions relative to their QED counterparts, 
particularly in the energy distribution, using these differences as a test of NCQED appears difficult because of 
their z^ and energy dependence. The search for ^-dependence remains the best possibility if the e"*" and e~ 
beams are unpolarized. 

There are, however, CP violating terms linear in z in the individual helicity amplitudes, as can be seen in 
Eq. itTTll or in Appendix [HI To probe these terms it is necessary to consider polarization effects. 
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3.2. Polarized Cross Sections 


We will confine our discussion to cross sections arising from longitudinally polarized beams. In this case, a 
typical cross section can be written 

^ [(1 + ^’e-)(l + -Pe+) + (1 - Pe-)(1 - Pe+) O-LL 

+ il + Pe-)a-Pe+)<yRL + il-Pe-)a + Pe+)<rLR] , (21) 

where, for example, denotes the cross section when the e~ beam has pure right-handed polarization {Pg- = 1 ) 
and the beam has pure left-handed polarization {Pe+ = —1)- The remaining cross sections are defined 
similarly. 

For the process e+e“ ^ 777 ) amplitudes with A = A vanish, and we can express the polarized cross section as 

ua 




{1-Pg-Pg+) 


CTRL + O-LR 

4 


Pe- - Pe+ f^LR - CTRL 
1 — Pe-Pe+ f^LR + CTRL 


(1 Fg-Pg+)(J,^npol [1 Pe¥i-^LR\ ) 


where the effective polarization PgS and the left-right asymmetry Arr are 


Peff 


Alr 


Pe- — Pe+ 

1 -Pg_Pg+ ’ 

C^LR — CTRL 
C^LR + CTRL 


( 22 ) 


(23) 

(24) 


The left-right asymmetry can be obtained using the squared amplitudes in Appendix o and Eq. which 
results in 


Alr = — 


El cos A 


(4C(3) 


29 47r2' 




AE"^ 

(log(^)-l)2(log(u;)-l) + 3 


(25) 


For the process e+e“ —> 777 , the NCQED correction is the main source of a left-right asymmetry. Competing 
standard model sources of left-right asymmetry such as Z exchange in Moller scattering are suppressed 
because they involve loops. Taking A^vc = 1.0 TeV, Tabled shows the cross section values for the cases 0 
Pg- = —Pg+ = ±0.9 and several values of and cos A = 1. 


yi TeV 

Co. 9 - 0.9 fb 

C-O.90.9 fb 

1.0 

4315.4 

4316.9 

2.0 

1187.9 

1189.3 

3.0 

557.1 

558.6 

4.0 

325.1 

326.6 

5.0 

213.9 

215.4 


TABLE I: The cross sections for 90% left-right and right-left polarized beams are shown for Anc = ITeV and A = 0. 

As the numbers in the Tabled indicate, the left-right asymmetry, though non-zero, must be distinguished from 
a fluctuation in the large left-right symmetric QED e+e" ^ 777 cross section. To obtain a sense of the range of 
values of Ajvc that can be probed polarized cross sections, we examined the signal to square root of background 
ratio 


R = 


P{o'0.9-0.9 — g'-O.QO.g) 
\/-^(o’o.9-o.9 + CT-o.go.g) 


(26) 
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FIG. 3: (color online) The bounds on Anc attainable using the left-right asymmetry of the total cross section are 
illustrated as a function of luminosity at y/S = 2 TeV (left) and as a function of y/S for £ = 500 fb~^ (right). The solid 
lines correspond to A = 0 and the dashed lines to A = n/4. 


Requiring R > 3 implies the bounds attainable on Ajvc illustrated in Fig.|21as function of luminosity C and 
collider energy y/S. 

The constraints on Anc obtainable from the polarized total cross section suggest that cuts on the polarization 
asymmetries in distributions such as 


{daLR- daRL)/dv {dcLR - daRc)/dcosO 

or 


(27) 


{daLR + daRL)/dv {da lr + da rl) /dcos 9 

could improve the bounds on Anc- These distributions are shown in Fig.^ While both distributions show a 


s 

b 



b 

+ 


s 

b 


n 


0.0 

COS0 


FIG. 4: (color online) The polarization asymmetries with respect to the photon energy fraction v (left) and the photon 
angle with respect to the beam axis (right) are shown. The three shaded regions correspond to center of mass energies of 
0.5, 1.0 and 5.0 TeV. 

distinct left-right asymmetry, the cos 9 distribution is the most promising from the experimental point of view in 
that it can be rather large - ^ few % - over a substantial region of cos 9. By imposing cuts on cos 9 it is possible 
substantially increase the lower bound on Anc obtained using Eq. The largest lower bound is obtained by 
restricting cos6> as | cos0| < 0.85, which is illustrated in Fig. 0. The behavior of the bound on Anc as the cut 
on cos 9 varies from | cos0| < 0.5 to | cos0| < 1.0 is shown in Fig.lHl 
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FIG. 5: (color online) Same as Fig.l^with |cos6| < 0.85. 
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FIG. 6: (color online) The bonnds on Ajvc for 0.5 < | cos^maxl < 1-0 are shown for y/S — 0.5 TeV (dot-dash-dash), 1.0 
TeV (dot-dash), 2.0 TeV (dashed) and 5.0 TeV (solid). Here £ = 500 fb”*^ and cos A = 1. The bounds scale as £^^^ and 
ycos A. 


4. DISCUSSION AND CONCLUSIONS 


To summarize, we have computed the noncommutative contributions to the angular and energy distributions 
of a single photon as well as to the total cross section for the process e“''e“ ^ 777 assuming both the unpolarized 
and polarized electron and positron beams. Because we are dealing with a three particle final state, it is possible 
to include these corrections using only the space-space portion of the tensor 0^''. This enables us to avoid the 
use of the space-time terms and thereby satisfy the requirement of unitarity 0 The use of space-time 
terms cannot be always avoided in 2 ^ 2 processes e.g. ey ^ ey, e+e” ^ 77 or e+e” ^ e“''e“ and this tends to 
complicate their interpretation. The cross sections and distributions depend on the angle A between the beam 
direction and the non-commutativity vector 9. In the unpolarized case, the noncommutative effects are second 
order in the ratio z = S/Aj^Q, whereas in polarized case, the noncommutative effects are leading order in z. 

In the unpolarized case, the shapes of the QED and NCQED energy distributions are quite different but it 
is the dependence of the cross section on the azimuthal angle 4> which offers the best opportunity to detect 
non-commutative effects. The observation of any variation of the cross section with respect to (/> is a clear 
violation of Lorentz symmetry. It is possible to introduce reasonable cuts significantly enhance this signature of 
non-commutativity. 

Further, the use of polarized beams makes it possible to probe the order z CP violating terms in the helicity 
amplitudes by measuring the left-right asymmetry. In contrast to Moller scattering, where Z exchange introduces 
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a large standard model left-right asymmetry which competes with the NCQED asymmetry, the NCQED left- 
right asymmetry in the 87 final state is the dominant source of asymmetry, with standard model contributions 
being suppressed by loops. Even without cuts on the polarized cross section, the bounds attainable on Ajvc 
are competitive with those obtained in pair annihilation |^. Imposition of cuts on cos0, the angle between one 
of the photons and the beam direction, extends the reach on Ajvc to the TeV range. Accumulating enough 
data to reach these bounds will require monitoring the (unknown) direction of the non-commutativity vector 
6. Techniques for doing this were proposed by Hewett, Petriello and Rizzo 0 and implemented by the OPAL 
collaboration 0 

Currently, the experimental lower bound on A^c is 140 GeV M , and the calculations of Ref indicate that 
Anc scales of 1.7 TeV can be probed in Moller scattering at a 500 Gev e+e” collider. Like Moller scattering, 
the NCQED contribution to e“'"e“ ^ 777 can be parameterized solely in terms of the unitarity preserving space- 
space components of 6^''. This, together with its NCQED dominant left-right asymmetry signature, makes the 
three photon process promising candidate in the experimental search for noncommutative effects. 
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APPENDIX A: SPIN AVERAGED CROSS SECTION 


The tensor 9ij can be parameterized in terms of a unit vector 9 and a noncommutativity scale A^rc as 

“ T2 ^ijk9 ■ 

^'■NC 


(Al) 


To define the coordinates, we fix the origin at the center of mass, choose the z axis parallel to pi and take 9 in 
the plane x-z plane. In this system, ki is defined by its polar angle ai, its azimuthal the angle /3i and its energy 
El. Similarly k 2 is defined by its energy E 2 , its polar angle a 2 and, for convenience, an azimuthal angle /3i - 1 -/ 32 ■ 

The phase space integration is given in detail in Appendix m where it is shown that, in addition to the 
variables mentioned above, it is necessary to introduce a minimum photon energy Emin to control the infrared 
singularities. Introducing the dimensionless variables {i = 1,2) 


El 

E 


Ci = cos ai 


En 


e = 


E 


n=\ I-^ . 


E 2 


the terms in Eq. m can be expressed as 


(A2) 


p = jvi (1 - nci) q = jU 2 (1 - nc 2 ) r = f (2 - (1 - nci) - 1 / 2 ( 1 - nc 2 )) 
s = jVi (1 -I- nci) t = ^V 2 (1 -I- nc 2 ) u= ^(2 — vi(\ + nci) — :/2 (1 + nc 2 )) 


(A3) 


The total cross section is expressible in terms of these variables as 


1 


(Jg+g— — 


6 (47r)^ 
/•i 


n2 — vi —e 


di/i 


di/2 


di/i 


di/i 


'l-vi 


L p2TT 

dci / dc 2 / I 
1 Jc- Jo 


dpi 




+7771 


V^(C+-C2)(C2-C_) 


X 


(A4) 
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where 

Cl (^1 +V2-1)T \/(1 - cf) (1 - i^i) (1 - 1 ^ 2 ) - 1) 

C± = Cl — z - 

1211/2 


(A5) 


Using the symmetry of Eq. under a permutation of ki,k 2 ,k 3 , and neglecting in the numerator of 
Eq. 11211 . the commutative contribution to the cross section becomes 


^QEd 

^ e+ e“ — >777 


2a3 

5 TT 


r* 1 — e ^2 — —€ 

di^i / d:>'2 + / di^i / dr's 

; 1 — J1 —e J € 


r‘=+ 


dci 


dC 2 


(2 — — ^ 2 )^ + {civi + € 21 / 2 )^ 

(1 - n^cl) (1 - n'^cl) ’ 


(A6) 


where the integration over (3i has been completed since the process has axial symmetry for 0 = 0. The integration 
of Eo. lAhll . neglecting terms which vanish for rr? jS 0, gives M 


jQed 

e+ g— —>.777 


2a3 


{\ogp- 1)^ (logw - 1) + (logp- l)logw + 3- ^ 


with 


S' 


P = 


U! = 


E^. 

mm 


(A7) 


(AS) 


As will be seen below, the NCQED correction is small relative to the pure QED cross section. To be certain 
that the comparison of the two is sensible, we computed the correction to the QED cross section obtaine d by 
retaining all the order terms in the numerators of Eq. eg. The calculations are explained in pdl IlfiL ll^ . 
Our result, which agrees with that of Berends and Kleiss, is 

2^3 


QED 


((logp-l)^(loga;-l) + 3) 


(A9) 


For the noncommutative term, the integrand is no longer invariant with respect to rotations about the z axis 
and it is necessary to consider the /3i integration in more detail. From Eq. EH we have 


ip = ^ki-0-k2 = \ {ki X ^2^ -0 


E 1 E 2 


— 2 — (sin ai (sin /3i cos Q !2 sin A + sin 122 sin /32 cos A) — cos ai sin Q !2 sin (/3i +132) sin A) , (AlO) 
2Anc 


where A is the angle between p) and 6. Assuming ip to be small, the noncommutative factor in Eq. d 
becomes = sin^(p = ip"^. Since /3i appears only in we can integrate this factor to obtain 


f -KZ 

/ d/3i(p^ = —[{l-‘icos^\){AciC2{vi+V2-l) + viV2{ci-C2)‘^) 

+8 cos^ A(1 - vi){l - V2){vi +V2- 1)] . 


(All) 


Owing to the additional factors of Ei and E 2 in Eq. ifTmli . the noncommutative contribution is infrared finite, 
and it is possible to set e = 0 in Eq. (IA4II . Then, using Eq. (Tmil and the symmetry of Eq. m in fci, ^ 2 , k^, the 
expression for the noncommutative contribution to the cross section is 


NCQED 
(J ^ — 

g-l-g 


64 TT S' 





M{ci,C2,Vi,V2) 
^C2Zf-, -T , 


(A12) 


where 


M{ci,C 2 ,Vi,l> 2 ) = [{2-Vi-V 2 )'^+n'^{ciVi+ 02 ^ 2 )"^] [{3, {vi + 122 ) - - V 1 V 2 ) - vlvl 

+ ( 2 i^iz ^2 ~ 7 )(i^i + 1/2)^ + 4 (j^i + 1/2)^ ~ iyi + ^2)'^ + {ciC2i2\V2(yi + 121122 3-122 — V) 

+ (1 — E2 + 12-^ + 121122 + 1^2 — ^j^i) + ^2^2 (7 — r'l 4 “ r'x 4 ” ^1^2 + ^2 — 2 r' 2 ))] (A13) 

[(1 - 3cl)l2il22 ( 4 ciC 2 (l^l 4 - 1^2 - 1) 4 - ^11^2(01 - C2)^) 4 - 8 c|(l - I 2 i)(l - V2){l2i +122- 1)] , 


X 
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and 


'D{ci,C 2 ,i'i,i' 2 ) = vlvl (1 - vi) (1 - V 2 ) {vi + - 1) (1 - ri^cl) (l - n^c|) sj (c+ - C 2 ) (c 2 - c_). (A14) 


When Ea. llA12ll is integrated over ai and 02 , we can take p = rr? jS = 0 except for terms involving logp, and 
take n = 1 in Eq. which does not introduce any collinear divergences. The integration of Eq. can be 

evaluated analytically and the noncommutative contribution is 


NcQED 
e+e~ —>777 


[2231 



^C(3)+sin2A(^^ + ^C(3) 


148957 71ogp\ 
34560 ~ 960 ) 


(A15) 


Erom this expression it is clear that production of the three photons depends on the angle between the incident 
beam and the vector 0, which violates Lorentz invariance. Since the direction of 9 is not known, we average over 
sin^ A to obtain 


-NCQED ^ f 65219 _ Tf__ 3 , X _ 7 log/3 \ 
o'e+e-^777 s V69120 480 4^^ ’ 1920 ) 

Adding Eq. and Eq. ITmli we get the total cross section 


^e+e“ 


^777 



2 (logp- 1)^ (logw 


l) + & + z^ 


/65219 
V69120 



7 log/3 \ 
1920 ) 


(A16) 


(A17) 


APPENDIX B: PHASE SPACE 


Taking the center of mass energies of the electron and positron to be E and the flux factor \v{ — V 2 \ ^ “I, the 
expression for the three body phase space in the electron-positron center of mass is 


dr = 


(27r) 


3! 2 (4£;2) 

1 


+P 2 -ki-k 2 - fcs) 


d^ki 


d^k 2 


d^k. 


35 (47r)® 


S{2E - El-E 2 -E 3 ) 5^ (/Ti -h k 2 + h'j 


{2Try 2Ei {2Try 2E2 {2 eY 2E3 
d^kid^k 2 d^k 3 


El E 2 E 3 


(Bl) 


where we have already introduced the symmetry factor 1/3!. In spherical coordinates the components of ki and 
k 2 can be defined by 


ki = El {sinai cos/3i, sinai sin/3i, cosai} 

fe = E 2 {sina 2 cos (/3i-I-/32),sinQ;2 sin(/3i-I-/32),cosa2} . 


(B2) 


The integration over <5^ (jsi + k 2 + ks^ gives 


\l + *2) 


E3 = y ks = 

= ijEf + E 2 + 2 E 1 E 2 (cos tti cos 02 + cos (^2 sin ai sin 02 ) 
where we have used the definitions (IB2II . The phase space integral Eq. m can then be written 

6i2E -E 1 -E 2 - Eg) d^kid^h 

E 1 E 2 


(B3) 


dr = 


35 ( 47 r )5 E1E2E3 

6{2E- El -E 2 - Eg) 


35 (47r)^ 


dEidE 2 d cos ogd cos a2d/3id/32 - 


Eg 


(B4) 
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with E 3 given by Ea. llB3ll . The remaining integration over 5{2E — Ei — E 2 — E^) leads to 

Eg = E 1 + E 2 2 Ei E 2 (cos «! cos 02 + cos (32 sin Oi sin 02 ) 

= { 2 E-Ei-E 2 f . 

Solving Eg. (IB 511 for (32 results in two solutions 

(32 = ±/320 

with 

/ EiE 2(1 — cosoi cosa 2 ) — 2E(Ei + E 2 — E)' 

P 20 = arccos --—^^- 

\ E 1 E 2 sinoi sina 2 

The derivative of the argument of the S with respect to ( 32 , 

d {2E — El — E 2 — E 3 ) E 1 E 2 . . . „ 

— ■ sinoi sina2 sinp2 . 


d(32 


E^ 


then gives 


(B5) 

(B6) 


(B7) 


(B8) 


S{2E -E 1 -E 2 


Es) = 


1 


E\E 2 „ 

——— sin ai sin a 2 sin P 20 
^3 


5{(32 — P20 ) + S {(32 + ( 32 o ) 


Replacing the Eq. fB^ in Eq. fB4ll we get 


(B9) 


dr 


1 

35 (47r)^ 


d£'idif 2 d cos Old cos a2<3(3iA(32 


^{(32 — (32o) + 5{(32 + (32o) 
sin oi sin 02 sin /32o 


In the squared amplitude, (32 is present only in the noncommutative factor . As can be seen from Eq. iTOl . 
the other variables are independent of (32- The integration over (32 produces w^(/32o) and w^{—(32o), but these 
contributions are equal after the subsequent (3i integration. Hence, for both the QED and NCQED terms, we 
can write 


dr 


2 dEli dif 2 dcosoi dcosa 2 d/3i 
(47r)^ sinoi sin 02 sin/32 


(BIO) 


where for simplicity we have dropped the zero on ^ 2 - The limits of integration on, say, cos 02 are constrained by 
Eq. IIB7II . Using this equation to solve for sinai sina 2 sin/32, we find, in the notation of Eq. (IA2II . 

sinoi sina 2 sin/32 = \/ (c+ — C 2 )(c 2 — c_) (Bll) 


with 


Cl (^1 + J^2 - 1) T \/(l - cj) (1 - vi) (1 - 1 ^ 2 ) (j^i +1^2-3) 

C± = Cl — z - 

UIU 2 


(B12) 


The limits of integration on the energies are determined by solving the relations 

33111111 ^ Ej < E El + E 2 + E^ = 2E , (B13) 

with j = 1, 2, 3, for Ei and E 2 . Elimination of E 3 yields the additional inequalities 

E 2 > E — El, < 2E — El — Emin • (B14) 

The limits of integration on E 2 depend on whether Emin < Ei < E — Emin or E — Emin < Ei < E. In the former 
case, 

E — El < E 2 < E , (B15) 

while in the latter 


Emin < E2 < 2E — El — Emin • 


(B16) 
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APPENDIX C: SQUARED MODULUS OF THE HELICITY AMPLITUDES 


The helicity amplitudes with permuted photon helicities can be derived by the corresponding permutation 
of the variables p, q, r and s, t, u. Amplitudes with every helicity reversed just change the sign of CP-breaking 
term, while amplitudes in which the electron and positron are exchanged can be derived changing the sign of 
the antisymmetric term and exchanging the variables p,q,r with s,t,u. The amplitudes {±, ±; Ai, A 2 , A 3 } and 
{A, A; ±, ±, ±1 are zero. The twelve non-zero helicity amplitudes are given in the following table. 


A, A; Ai, A 2 , As 

1-^^ A, A; Ai, A 2 , A 3 1 

A, A; Ai, A 2 , As 

1-'^^ A, A; Ai, A 2 , A 3 1 

+i +! +! “ 

pqst 

+ 

J. 

J 

J 

+ 

^A- 

pqst 

+ 

J. 

+ 

J 

+ 

A 

prsu 

+ 

J. 

J 

+ 

1 

r 

A- 

prsu 

+1 +) + 

J^A- 

qrtu 

+ 

J. 

+ 

J 

1 

qrtu 

J 

+ 

J 

J 

+ 

pqst 

—, -f; -f, -f, — 

pqst ^ 

J 

+ 

J 

+ 

1 

—A+ 

prsu 

“1 +; +! “! + 

. 

prsu 

+; +! “! “ 

-^A+ 

qrtu 

“1 +; “! +! + 

— 

qrtu 


TABLE II: The absolute squares of the twelve non-zero helicity amplitudes, |A1 a,A;Ai,A 2,A3 l^> expressed in terms of 
the amplitudes A± given below. 

The common factor A± in all the entries in Tabled is 


A± = 2 e® 


5 - 4^2 -S' + 


a{p'^ + s^) b{q^ + t^) c{A+u'^) 


2 ' be ac ab 

pq+ {p + s){q + t) + st pr + {p + s){r + u) + su qr + {q + t){r + u) + tu 


± ^vw ( - + 7 + - ) £{ki,k2,Pi,P2)) 

'a b c 


(Cl) 
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